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SUMMARY 


Extensive wind tunnel tests were conducted to establish the preflight aero- 
dynamics of the Shuttle vehicle. The venicle was designed and built using this 
data. In an effort to verify the preflight aerodynamics, a flight test program was 
established. Maneuvers were made and data taken at various angles-of-attack and 
Mach numbers throughout the descent flight envelope. 

This paper presents the longitudinal short period aerodynamics of the space 
shuttle Columbia as determined from flight test data. These flight-determined 
results were compared with the preflight predictions. In spite of the scatter in 
the flight-determined parameters, the general trends of the parameters with Mach 
number led to the conclusion that with the exception of the control effectiveness in 
the vicinity of Mach=l, the preflight predictions were representative of the 
vehicle's short-period aerodynamics. The pitch RCS was found to be more effective 
than predicted. 
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INTRODUCTION 


The Space Shuttle vehicle has received one of the most extensive preflight 
analysis of any aircraft that has ever flown. Thousands of wind-tunnel hours went 
into its development and refinement. The description of the Shuttle aerodynamics 
that resulted from the wind-tunnel tests and the analytical studies is a detailed 
description of the aerodynamic characteristics of the Shuttle vehicle over a flight 
envelope covering a Mach number range from 27 to 0 (reference 1). In an effort to 
verify the preflight aerodynamics, a flight test program was established. The 
program was planned as an on-going process based on the analysis of measurement data 
from each succeeding flight. Since only a limited number of maneuvers could be 
performed during a given Shuttle descent, these were planned to examine different 
aspects of the Shuttle aerodynamics so that as much of the flight envelope as 
possible could be verified. 

The difficulty with this plan-of-attack was that the types of maneuvers that 
could be performed, within the constraints of safety and that were allowed by the 
Shuttle flight control system, were not the best maneuvers for identifying the 
Shuttle aerodynamic parameters. However, since the maneuvers that were allowable 
represent the only data available for analysis, each data set was examined in 
extreme detail. The majority of the results shown herein were obtained by 
processing the test maneuvers with a Maximum Likelihood parameter extraction 
program. Also, where appropriate, selected runs were examined using analytical 
techniques. When these multiple analyses show similar results, the confidence in 
the values determined Is Increased. 

This paper will represent the results of analyzing the longitudinal maneuvers 
from five Columbia flights. These results will be compared with those of 
reference 1, and the Implications and significance of any differences will be 
discussed. Also, comments will be made as to the confidence in the parameter values 


SYMBOLS 

acceleration measured along X, Y, and Z body axes, 
respectively, g units 

wing span, m (ft) 

wing mean geometric chord, m (ft) 

force along X, Y and Z body axes, respectively, N (lb) 

acceleration due to gravity, m/sec^ (ft/sec^) 

moment of inertia about X, Y, and Z body axes, 
respectively, kg-m^ (slug-ft^) 

product of inertia, kg-m2 (slug-ft^) 

% 

distance from airplane center of gravity to center of pressure 
of horizontal tail, m (ft) 

rolling, pitching, and yawing moments, respectively, N-m 
(ft-lb) 
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mass, kg (slugs) 

rate of roll, pitch and yaw, rad/sec or deg/sec 
dynamic pressure, N/m^ (slug/ft^) 
estimate of error covariance matrix 
wing area, m^ (ft^) 

velocity along X, Y, and 7 body axes, respectively, m/sec 
(ft/sec) 

control vector 

airplane total velocity, m/sec (ft/sec) 

matrix of measured states and input variables 

body coordinate axes through airplane center of gravity 

state vector 

output vector 

me a s u r ement vector 

angle of attack, rad or deg. 

angle of sideslip, rad or deg. 

elevator deflection, rad or deg. 

speedbrake deflection, rad or deg. 

pitch angle, roll angle, rad or deg. 

parameter vector 

air density, kg/m^ (slug/ft^) 

measurement noise vector 

perturbation in parameter vector 

pitching-moment coefficient, M^/ qSc 


axial-force coefficient, F /qS 
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Subscripts : 

computed 
index 
measured 

coefficient at trimmed conditions 
trimmed cond i t ions 

forces in moments due to Reaction Control System (RCS) 
moments due to aerodynamic parameters 
inertial moment 

Superscripts : 

-1 inverse matrix 

T transpose matrix 

M measured quantity 

0 nominal evaluation 

estimated value 

A dot over a symbol signifies a derivative with respect to time* 
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TEST VEHICLE 

The orbiter configuration is shown as figure 1, and key physical characteris- 
tics are given in Table I • The thick, double delta wing is configured with full 
span elevons, comprising two panels per side* Each elevon panel is independently 
actuated* All four panels are deflected symmetrically as an elevator (6 e ) for 

pitch control and left and right elevons are deflected differentially as an aileron 

(6 a ) for roll control* 

The body flap is used as the primary longitudinal trim device* The elevons are 

programmed to follow a set schedule with the body flap deflection to provide a 

desired aileron effectiveness. 

The vertical tail consists of the fin and split rudder. The rudder panels are 
deflected together for yaw control and are separated to act as a speedbrake (6 SB ) 
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to provide for subsonic energy modulation. The speedbrake opens fully (87.2 
degrees) just below Mach 10 and then follows a predetermined schedule until Mach 0.9 
is reached. The rudder is not activated until until Mach 3.5. 

Stability augmentation is provided by the aft reaction control system (RCS) 
jets, with the forward jets reserved for on-orbit attitude control and for aborts. 
The aft yaw jets are active until Mach 1, while the pitch and roll jets are termi- 
nated at dynamic pressures of 20 and 10 psf, respectively. Additional details on 
the Shuttle vehicle and its systems are given in reference 1. 


MANEUVERS 

During STS-2 through 5 and STS-9, specially designed maneuvers were performed 
to obtain data for use in the parameter extraction programs. These maneuvers were 
performed to obtain data at specific points during the descent trajectory. The test 
points were chosen so that aerodynamic parameters could be determined along the 
descent trajectory to verify the aerodynamic model obtained from the wind tunnel 
tests. This verification procedure will add confidence to the assumed aerodynamics 
of the Shuttle where there is agreement and will point to areas of potential 
inaccuracy where there is no agreement. 

The actual forms of the inputs to be performed were developed using a Shuttle 
simulation to generate responses for various inputs and then extracting parameters 
from these responses. The control inputs that gave the best definition of the 
parameters of interest were then used for the flight tests. In spite of the care 
taken to design effective inputs, since the automatic control system was active, the 
controls were coupled and the resulting responses were reduced in magnitude and 
correlated with each other and the control inputs. This led to ident if lability 
problems and correlation of parameters during the extraction process. Additional 
details on the maneuver design are given in reference 2. 


INSTRUMENTATION AND DATA PROCESSING 

As a development vehicle, the Shuttle is fully instrumented and has a number of 
redundant systems for measuring various vehicle states and controls. Three Instru- 
ment packages will be mentioned specifically. First is the Aerodynamic Coefficient 
Identification Package (ACIP), an instrumentation package specifically designed to 
measure rates, accelerations and control surface positions required for parameter 
identification. The ACIP data were recorded at 172 samples-per-second . Second is 
the instrumentation for the flight guidance and control system (RGA, AA) which Is a 
source for acceleration and rate measurements. The RGA, AA data were recorded at 25 
samples-per-second, but are very noisy. The third source of flight measurements is 
the navigation Instrumentation (IMU). The IMU measurements are high fidelity, but 
are only recorded at one sample-per-second , which limited their usefulness. 

With the exception of STS-2 , where data were not available because of a 
recorder failure, the ACIP data were the primary source for the linear and angular 
accelerations, angular rates, and control surface deflections. The RCS chamber 
pressures were used to determine the jet thrust, and these measurements came from 
the vehicle operational instrumentation. 

The data considered most reliable were used to generate a best estimated 
trajectory (BET) for the Shuttle vehicle. The data put on the tapes prepared for 
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parameter extraction consisted of only those maneuvers considered appropriate for 
extraction. The linear and angular accelerations, angular rates and control surface 
deflections came from the ACIP instrumentation except for STS-2 where a combination 
of IMU and RGA, AA data was used. The BET angular rates and linear accelerations at 
the start of a maneuver were taken as initial conditions, and the rates and 
accelerations were integrated over time to obtain angular positions and vehicle 
velocities. The velocities were then corrected for the effect of winds and the 
resulting components were used to calculate the vehicle total velocity, angle- 
of-attack and angle-of-sideslip. This combined data set comprises the data 
contained on the tape to be processed by the parameter extraction program, and is 
recorded at 25 samples-per-second. Additional details on the instrumentation and 
data processing can be found in references 3, 4, and 5. 


EXTRACTION METHODS 

Several methods were used to estimate stability and control parameters using 
the flight test data. These were a maximum likelihood technique described in 
reference 6, a regression technique described in reference 7, and two analytical 
techniques described in references 8 and 9. The Maximum Likelihood technique 
utilizes the log-likelihood function 

j(6) = - J l R 1 ri ± - -| log |r| 


where 


- 3y i 

"i - Z 1 - h ■ y i - y i ( V ' W 


A0 

e o 


with the measurement vector and y^ the output vector which comes from 
x«f(x, U, 0 , t) and y=g(x, U, 0, t) . In the above equation is assumed to have 
a Gaussian distribution and the representation x=f(x, U, 0, t) is assumed to 
accurately represent the physical system. The unknowns to be estimated are the 
elements of 0 and R. Minimizing J with respect to R, 
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is obtained. The estimates for the parameters are obtained from the equation 
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yielding the parameter estimates 9=9 q+A0. 

The regression technique utilizes the cost function 

N 

J (0 ) = y [ x . - f .(x, U. 9 ) ] 
r ■- ri ri ' r ' 
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where r indicates the rth state equation. The estimates of the unknown parameters 


are obtained from the equation 
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0 which results in 
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where matrix x^ includes measured states and output variables (assumed noise 
free). Another feature of the stepwise regression program is that it estimates the 
parameters according to the percent of the vehicle motion explained by that partic- 
ular parameter, identifying the most important parameters in the mathematical model. 
The parameter importance feature has worked reasonably well for the Shuttle data, 
but because of the reduced signal-to-noise ratio of the flight data, the parameters 
extracted using the regression program were biased in many cases. 


Two analytical techniques were used to supplement the analysis of the 
data using the extraction computer programs. One analytical technique was 
determine the longitudinal static parameters. Assuming that angular rates 
constant , 
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(where C^, C L and the states are determined for a particular angle-of-attack 
(a)) were used to calculate and for specific angle-of-attack values. Then 
as angle-of-attack varied, so did and These variations were used to 

calculate AC m /Aa and AC^/Aa. 


Another analytical technique used was deterministic analysis. The form of the 
equation used is: 


Vehicle total moment = RCS moment + aerodynamic moment + unmodeled moment. 

The deterministic equations can be written for moments about any of the vehicle body 
axes. The relative magnitudes of the terms that make up the aerodynamic moments 
indicated which of the aerodynamic parameters have the most influence on the vehicle 
total moment. The extent of the match of shape and magnitude between the aero- 
dynamic moment and the vehicle total moment is an indication of how well the assumed 
mathematical model will represent the measured vehicle motion. 


EQUATIONS OF MOTION 

The equations used in this program are perturbation equations from trimmed 
level flight and are written relative to the set of body axes shown in figure 1. 

The equations used to describe the longitudinal motions were 
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RESULTS AND DISCUSSION 

The allowable inputs for parameter extraction maneuvers were small and so the 
data were marginal for use with the parameter extraction programs. Because of this, 
several methods of analysis were used in the hope that there would be agreement in 
the parameter values determined from several sources and this would add confidence 
to the results. The stepwise regression program was also tried as an analysis tool, 
but the response to input was so small that results using the regression program 
were not usable. Therefore, the analysis procedures actually used to examine the 
flight data were a Maximum Likelihood program and two analytical techniques. 

Only a limited number of the longitudinal maneuvers were usable for parameter 
extraction. The push-over, pull-up maneuver or the pull-up, push-over maneuver 
proved to give the most consistant results and yield parameter values with the 

lowest estimated standard deviations. Even with this "best" maneuver, only C 7 , 

/j ct 

C^, and Cmgg were identifiable. C Z( j e was not considered identifiable 
since the extracted values exhibited considerable scatter, and the estimated 
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standard deviations were over 30 percent of the estimated parameter value* 
Therefore, was at reference l values, except in the vicinity of 

Mach^l , for the runs shown in this paper* The parameter values chosen for 

are given in the sections discussing the individual parameters* 



The values determined for , and are given in Table II 

for various Mach numbers* The estimated standard deviations are also given In the 
table* These values are an Indication of confidence in the extracted parameter 
value. An estimated standard deviation that Is less than one-tenth of the value of 
the estimated parameter indicates some confidence that the estimated parameter is 
reasonable. The parameters values for C 7 , C m , (L r , and C 7 * are 

plotted versus Mach number. There are two plots on each figure, one for Mach > 5 
and one for Mach < 5. The individual parameters will now be discussed. 


'%<x 


- The variation of normal force coefficient with angle-of-attack 


parameter is shown plotted versus Mach number as figure 2. Examination of the 
figure Indicated a trend toward values that were more negative than those of 

reference 1 for Mach 10 and above. Below Mach 10, extracted values of C 7 tend 

L a 

to be less negative than those of reference 1. In this lower Mach range, several 

parameter values are shown at the same Mach numbers. Because C 7 p was not con- 

^<5e 

sldered identifiable from the runs available, selected values were used for this 
parameter In the vehicle mathematical model. Using reference 1 values for Mach 
numbers above 2 resulted in reasonable values for the other parameters. However, 
between Mach 1 and Mach 2, reference 1 indicated a change in the value of C 7 

by a factor of 5. When these larger values were used in the extraction mathematical 
model, the values for the other extracted parameters showed greater scatter and in 


some cases had unreasonable magnitudes • Because the reference 1 values for C z 

resulted in unreasonable values for CL and Cr 

™a 1 

per radian to .6 per radian range were chosen (see figure 5). 
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The values calculated for C z ^ , using analytical techniques and assuming 

linear variations in alpha, "Z"-body acceleration and elevon deflection and pitch 
rate constant are also shown on figure 2. These values generally showed trends 
similar to those of reference 1. Below Mach 7 the magnitude of the values 
calculated for tended to be less negative than those of reference 1. In 

most cases, the calculated values showed the same trends as the values determined 
using the maximum likelihood extraction program. 


^ma ~ stat * c stability parameters for the flights discussed in this 

paper, plotted versus Mach number, are shown as figure 3. The values from 
reference 1 are designated by a solid line. Especially at the higher Mach numbers, 
there are many regions where changes by 40 to 50 percent in a very small 

range of Mach numbers. These abrupt variations make comparing trends and values 
between the extracted values and those of reference 1 difficult. In 

general, the extracted parameters seemed to show less static stability above Mach 10 
and more below Mach 10 than the preflight predictions. 


8 



The values of calculated using analytical techniques for the regions of 

the trajectory where alpha, "Z"-body accelerations and elevon deflection changes 
were linear and pitch rate was constant showed trends and magnitudes of the 
parameter values that were similar to those of the parameters extracted using the 
Maximum Likelihood program. This agreement tended to give additional confidence to 
the validity of the values determined. 

Cjjj^e — The control effectiveness parameter is plotted versus Mach number 

and shown as figure 4. The values determined using the Maximum Likelihood extrac- 
tion program generally showed the same trend but slightly less effectiveness than 
reference 1 down to Mach 5. Below Mach 5 the extracted parameters showed similar 
trends, but some of the values determined were different. The predominant feature 
in the preflight longitudinal aerodynamics is the abrupt change in the parameter 

values between Mach 1.5 and .5. In the case of CL. the effectiveness is 

^6e 

predicted to increase by a factor of 3 in the vicinity of Mach 1. Because of lack 
of maneuvers that would yield identifiable parameters during this part of the 
flight, this variation In effectiveness could not be verified. The trends of the 
extracted parameter values seemed to imply an increased effectiveness, but not as 
great an increase as predicted by reference 1. 

Cz^ e * The variation of force along the M Z"-body axis with elevon deflec- 
tion parameter from reference 1 and fixed values used for extraction of other 
parameters are plotted versus Mach number in figure 5. Above Mach 2 the values from 
reference 1 were used for However, below Mach 2 the plotted values seemed 

to give more realistic values for C^, and 

At Mach- .7 the results from the maneuvers were not as sensitive to variations 
in as those of Mach-. 6 and Mach-. 8. When Cz^ was set at the value 

from reference 1, became significantly positive, but had a very 

reasonable value. When was set at *^.4 per radian C^, and 

Cm< 5 e did not va ty significantly, however, the values extracted for 

showed much less control effectiveness than was expected. 

The pitch jets on the Shuttle were only active for dynamic pressures less than 
20 psf. The evaluation of the effectiveness of the pitch jets was made using 
several methods so that the values obtained could be verified * Values for the pitch 
jet effectiveness were obtained at three Mach numbers or alternatively at three 
values of q. Values were obtained using the Maximum Likelihood extraction program 
and deterministically from the equation 

* l y - < r Z - I X> ? r - (r2 - » 2) r XZ * M y + \ P V 2 S? [ C ° + C m l 1 

a 6e 

where C- and L. were obtained from reference 1. The results of this 
ra a m 6 e 

determination of jet effectiveness is shown as figure 6. The figure shows the 
percent effectiveness plotted against dynamic pressure. The two methods agreed 
within 10 percent. Indicating a strong confidence in the results. The flight- 
derived values, however, are very different from the predicted values at the lower 
values of dynamic pressure but tend to agree better with the predicted values at the 
higher values of dynamic pressure. 
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The results show that the up-firing jets are about 15 percent more effective 
than the down-firing jets. This is a measure of pitch jet interactions, and it had 
been hoped that by working in the dynamic pressure regime near zero, the interac- 
tions due only to impingement could be determined. Then by comparing the results to 
pitch jet effectiveness and higher dynamic pressures, the flow field interactions 
could be separated from the impingement interactions. Results from this attempt 
were inconclusive because of instrumentation resolution at low dynamic pressures and 
by the fact that at 20 psf, when jets were cut off, the flow effects were poorly 
defined . 


CONCLUDING REMARKS 

Longitudinal aerodynamic parameters were extracted from flight test data for 
five flights of the Space Shuttle Columbia. These parameter values were compared 
with the preflight predictions to try to determine how well the predictions 
described the Shuttle aerodynamics. In general, the values of the identifiable 
parameters exhibited trends with Mach numbers that were similar to those of 
reference 1. However, for most of the Mach range, the actual values were different. 

The extracted values of C£ a were more negative than those of reference 1 

for Mach numbers greater than 10 and less negative for Mach numbers less than 10. 
However, the correlation problems that resulted from the small magnitudes of the 
vehicle’s response to inputs reduced the confidence in the low Mach number results. 

The extracted values for tended toward less static stability above 

Mach 5 than those of reference 1. Below Mach 5, the parameter values extracted from 
flight data agreed reasonably well with those of reference 1; however, lack of data 
prevented a good definition of the large variation In around Mach 1. 

The extracted values for seem to generally agree with those of 

reference 1 in the Mach 1.5 to Mach 20 range. In this range, the parameter was 
fairly well identified, lending some confidence to the values obtained. Above 
Mach 20, the parameter was not well identified, and considerable scatter was evident 
in the values determined. Below Mach 1.5, the values determined indicated less 
elevon effectiveness than predicted by reference 1. When the reference 1 values 
were used for c z<$ e t ^ e same Mach range, C% and had unreasonable 

values. If the reference 1 values of C 7c were cut in half, then the C 7 

and values were more reasonable. 

The RCS pitch jet effectiveness determined from flight indicates that the pitch 
jets are approximately twice as effective as predicted over most of the range where 
they are active. 

With the exception of control effectiveness in the vicinity of Mach=l, the 
majority of the flight conditions covered during the entry phase of the five 
Columbia flights analyzed, the agreement with the preflight predictions is close 
enough to conclude that these predictions are a reasonable representation of the 
vehicle’s short period aerodynamics. 
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TABLE I.- ENTRY PHYSICAL CHARACTERISTICS OF SPACE SHUTTLE ORBITER 

Mass propoerties (range for five flights): 


Mass, kg 91,917-100,309 

Moments of Inertia (range for five flights): 

l x , kg-m 2 1,171 ,428-1 ,313,633 

2 

l y , kg-m 9,228,939-9,614,705 

2 

l z , kg-m 9,584,958-10,031,878 

2 

! X z> 205,832-223,189 


Wing: 

2 

Reference area, m . . • 249.91 

Mean aerodynamic chord, m 12.06 

Span, m 23.79 


E levon (per side ) : 

2 

Reference area, m 19.51 

Mean aerodynamic chord, m . . . . 2.30 


Rudder (per side panel) : 

2 

Reference area, 9.30 

Mean aerodynamic chord, m ....... 1.86 


Body Flap: 


Reference area, 12.54 

Mean aerodynamic chord, m . ♦ ..... 2.06 
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Maximum Likelihood Maximum 
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TABLE II.- EXTRACTED PARAMETERS (CONTINUED) 
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TABLE II.- EXTRACTED PARAMETERS (CONTINUED) 
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Figure 2.- Change in normal force with alpha versus Mach number. 
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Figure 3.- Static stability versus Mach number. 
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Figure 5.- Change in normal force with elevon deflection 
versus Mach number. 
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RCS PITCH JET EFFECTIVENESS 
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Figure 6*- Pitch jet effectiveness versus dynamic pressure. 
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